Quantum oscillations in electron doped high temperature superconductors 
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Quantum oscillations in hole doped high temperature superconductors are difficult to understand within the 
prevailing views. An emerging idea is that of a putative normal ground state, which appears to be a Fermi liquid 
with a reconstructed Fermi surface. The oscillations are due to formation of Landau levels. Recently the same 
oscillations were found in the electron doped cuprate, Nd2-xCexCu04, in the optimal to overdoped regime. 
Although these electron doped non-stoichiometric materials are naturally more disordered, they strikingly com- 
plement the hole doped cuprates. Here we provide an explanation of these observations from the perspective of 
density waves using a powerful transfer matrix method to compute the conductance as a function of the magnetic 
field. 



I. INTRODUCTION 

Periodically new experiments tend to disturb the status 
quo of the prevailing views in the area of high tempera- 
ture cuprate superconductors. Recent quantum oscillation 
(QO) experiments' fall into this category.'' The first set 
of experiments were carried out in underdoped high quality 
crystals of well-ordered YBa2Cu306+5 (YBCO), stoichio- 
metric YBa2Cu408 (Y124) and the overdoped single layer 
TlaBaaCuOe+s"*. 

More recently oscillations are also observed in electron 
doped Nd2-xCexCu04 (NCCO)." The measurements in 
NCCO for 15%, 16%, and 17% doping" are spectacular The 
salient features are: (1) The experiments are performed in the 
range 30 — QAT, far above the upper critical field, which is 
about lOT or less; (2) the material involves single CuO plane, 
and therefore complications involving chains, bilayers, Ortho- 
II potentialr^ etc. are absent; (3) stripes'^' may not be germane 
in this case.'^ It is true, however, that neither spin density 
wave (SDW) nor d-density wave (DDW)'^ are yet directly ob- 
served in NCCO in the relevant doping range, but QOs seem 
to require their existence, at least the field induced variety (see, 
however Ref.[T6h; (4) these experiments are a tour de force be- 
cause the sample is non-stoichiometric with naturally greater 
intrinsic disorder. The effect is therefore no longer confined to 
a limited class of high quality single crystals; (5) The authors 
have also succeeded in seeing the transition from low to high 
frequency oscillations'^ in NCCO as a function of doping. 

Here we focus on NCCO. We shall see that disorder plays 
an important role. Without it it is impossible to understand 
why the slow oscillations damp out below SOT for 15% and 
16% doping, and below GOT for 17% doping, even though the 
field range is very high. For 17% doping, where a large hole 
pocket is observed coiTesponding to very fast oscillations (in- 
consistent with any kind of density wave order), the neces- 
sity of such high fields can have only one explanation, namely 
to achieve a sufficiently large uJcT, where ojc = eB/m*c, r 
is the scattering lifetime of the putative normal phase, m* 
the effective mass, and B the magnetic field. Qualitatively, 
the Dingle factor, D, that suppresses quantum oscillations is 
D = e~P^/'^<^'^ where p is the index for the harmonic. As- 
suming a Fermi velocity, suitably averaged over an orbit to be 
Vf, the mean free path / = vpT. Thus D can be rewritten 



as D 



-pTihckf / eBl 



A crude measure for fc^? is given by 



expressing the area of an extremal orbit. A, as A = Trkp. By 
setting m*VF — hkp the explicit dependence on the parame- 
ters m* and vp was eliminated. Assuming that the mean free 
paths for the hole and the electron pockets are more or less 
the same, not an unreasonable assumption, the larger pockets, 
with larger kp, will be strongly suppressed for the same value 
of the magnetic field because of the exponential sensitivity of 
D to the pocket size. This argument is consistent with our ex- 
act transfer matrix calculation using the Landauer formula for 
the conductance presented below. 

Here we show that the oscillation experiments in NCCO 
reflect a broken translational symmetry^^ that reconstructs the 
Fermi surface in terms of electron and hole pockets.^ The em- 
phasis is not the transfer matrix method itself, but its use in 
explaining a major experiment in some detail. We study both 
SDW and singlet DDW orders with the corresponding mean 
field Hamiltonians. A more refined calculation, beyond the 
scope of the present paper, will be necessary to see the subtle 
distinction between the two order parameters. 

In Sec. II we introduce our mean field Hamiltonians and in 
Sec. Ill we discuss the transfer matrix method for the com- 
putation of quantum oscillations of the conductance. Sec. IV 
contains the results of our numerical computations and Sec. V 
our conclusions. 



II. MEAN FIELD HAMILTONL4N 

We suggest that the experiments in NCCO can be under- 
stood from a suitable normal state because the applied mag- 
netic fields between 30-65 T are so far above the upper critical 
field, which is less than 10 T, that vortex physics and the su- 
perconducting gap are not important. Our assumption is that 
a broken translational symmetry state with an ordering vector 
Q = (tt/o, 7r/a) (a being the lattice spacing) can reconstruct 
the Fermi surface resulting in two hole pockets and one elec- 
tron pocket within the reduced Brillouin zone, bounded by 
the constraints on the wave vectors ± ky = in/a. One 
challenge here is to understand why the large electron pock- 
ets corresponding to 15 and 16% doping resulting from the 
band structure parameters for NCCO defined below are not 
observed, but the much smaller hole pockets are. Another 
challenge is to understand why the large Fermi surface at 17% 
doping is not observed until the applied field reaches about 60 
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T. The reason we believe is the existence strong cation dis- 
order in this material. It is therefore essential to incorporate 
disorder in our Hamiltonian. For the Hamiltonian itself, we 
consider a mean field approach, and for this purpose we con- 
sider two possible symmetries, one that corresponds to a sin- 
glet in the spin space (DDW) and one that is a triplet in the 
spin space (SDW). Note that these are particle-hole conden- 
sates for which orbital function does not constrain the spin 
wave function unlike a particle-particle condensate (supercon- 
ductor) because there are no exchange requirements between 
a particle and a hole. 

We believe that it is reasonable that as long as a system is 
deep inside a broken symmetry state, mean field theory and 
its associated elementary excitations should correctly capture 
the physics. The fluctuation effects will be important close to 
quantum phase transitions. However, there are no indications 
in the present experiments that fluctuations are important. The 
microscopic basis for singlet DDW Hamiltonian is discussed 
in some detail in Refs. 19.20 and in references therein. So, 
we do not see any particular need to duplicate this discussion 
here. The mean field Hamiltonian for the singlet DDW in 
real space, in terms of the site-based fermion annihilation and 
creation operators of spin a, Ci^a and cj^, is 



H 



DDW 



Ci,, 



+ h.c., (1) 



where the nearest neighbor hopping matrix elements are 

tu-.. = (2) 

iWn 

(3) 



Here Wq is the DDW gap. We also include the next nearest 
neighbor hopping t', whereas the third neighbor hopping t" is 
ignored to simplify computational complexity without losing 
the essential aspects of the problem. The parameters t and t' 
are chosen (see Table [III to closely approximate the more con- 
ventional band structure, as shown in Fig.[T] We have checked 
that the choice t" ~ provides reasonably consistent results 
for the frequencies in the absence of disorder For example, 
for DDW, and 15% doping, the hole pocket frequency is 185 
T, and the corresponding electron pocket frequency is 2394 T. 



Similarly, the SDW mean field Hamiltonian is 

HsDW = E + aVs{-iy-+'^] cl^cu. 



(4) 



and the spin cr = ±1, while the magnitude of the SDW am- 
plitude is Vs. In both cases a constant perpendicular mag- 
netic field B is included via the Peierls phase factor ai j — 

^ A • dl, where A = (0, —Bx, 0) is the vector potential 
in the Landau gauge. We note that usually a perpendicular 
magnetic field, even as large as 60T, has little effect on the 



TABLE I: The band parameters, the chemical potential, and the mean 
field parameters for DDW and SDW used in our calculation. F in 
Tesla corresponds to the calculated oscillation frequencies of the hole 
pocket, the so-called slow frequencies. The measured F for 15% 
doping is 290 ± 10 T and for 16% doping is 280 ± 15 T. The calcu- 
lated magnitude of F does depend on the neglected t" . 
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DDW 16% 
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SDW 15% 
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FIG. 1: (Color online) The solid curve represents the t — t' — t" 
band structure (t = 0.38eV, t' = 0.32t, t" = 0.5t'), and the 
dashed curve corresponds tot — t' band structure, (see Table |D. The 
quasiparticle energy is plotted in the Brillouin zone along the triangle 
(0,0) — s> (tTjO) — s> (7r,7r) — ^ (0,0). In the inset the chemical 
potential, fi, was adjusted to obtain approximately 15% doping. 



DDW gap,^' except close to the doping at which it collapses, 
where field induced order may be important. 

We have seen previouslyi^ that the effect of long-ranged 
correlated disorder is qualitatively similar to white noise inso- 
far as the QOs are concerned. The effect of the nature of disor- 
der on the spectral function of angle resolved photoemission 
spectroscopy (ARPES) was found to be far more important. 
The reason is that the coherence factors of the ARPES spec- 
tral function are sensitive to the nature of the disorder because 
they play a role similar to Wannier functions. In contrast, the 
QOs are damped by the Dingle factor, which is parametrized 
by a single lifetime and disorder enters in an averaged sense. 

Thus, it is sufficient to consider on-site disorder The on- 
site energy is (5-correlated white noise defined by the disorder 
average ei = and elej = Vq^i,}- For an explicit calculation 
we need to choose the band structure parameters, Wq, Vs, and 
the disorder magnitude Vq. When considering the magnitude 
of disorder one should keep in mind that the full band width 
is 8t. The magnetic field ranges roughly between SOT and 
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64T, representative of the experiments in NCCO. The mag- 
netic length is Is = y/h/eB, which for B = SOT is approxi- 
mately 12a, where the lattice constant a is equal to 3.95A. 

The effect of potential scattering that modulates charge den- 
sity is indirect on two-fold commensurate SDW or DDW or- 
der parameter, mainly because SDW is modulation of spin 
and DDW that of charge current. Thus, the robustness of these 
order parameters with respect to disorder protects the corre- 
sponding quasiparticle excitations insofar as quantum oscil- 
lations are concerned, as seen below in our exact numerical 
calculations. Thus we did not find it important to study this 
problem self consistently. 



proof of this curious fact. Clearly, larger eigenvalues con- 
tribute insignificantly to the more general formula for the 
conductance:^"^ 



2M 



a{B) = -TvY r 



^ ^(7ivr;) + (w;)-i + 2 



(8) 



When the eigenvalues do come in pairs, the conductance for- 
mula simplifies to the more common Landauer formula:— 



(9) 



III. TRANSFER MATRIX METHOD 

The transfer matrix method and the calculation of the Lya- 
punov sketched elsewhere'^ is fully described here for the 
case of singlet DDW; for SDW the generalization is straight- 
forward, where the diagonal term must be modified because 
of Vs, and the term Wq will be absent. Consider a quasi-lD 
system, L ^ M, with a periodic boundary condition along 
y-direction. Let = (i/'n,i' V'n,2, • ■ • , V'h.m)^ be the ampH- 
tudes on the slice n for an eigenstate with a given energy, then 
the amplitudes on three successive slices satisfy the relation 



Tn An T„ Bn 
1 



(5) 

where T„, A^, B^ are M x M matrices. The non-zero matrix 
elements of the matrix An are 



(An) m. 
[An) 

(-^n)m,m— 1 



4 
4 



(6) 



The transfer matrix method is a very powerful method and 
the results obtained are rigorous compared to ad hoc broaden- 
ing of the Landau levels, which also require more adjustable 
parameters to explain the experiments. Once the distribution 
of disorder is specified there are no further approximations. 
We note that the values of M were chosen to be much larger 
than our previous work,'^ at least 128 (that is 128 a in phys- 
ical units) and sometimes as large as 512. The length of the 
strip L is varied between 10^ and 10®. This easily led to an 
accuracy better than 5% for the smallest Lyapunov exponent, 
7;, in all cases. 

We have calculated the a6-plane conductance, but the mea- 
sured c-axis resistance, R^, is precisely related to it, at least 
as far as the oscillatory part is concerned. This can be seen 
from the arguments in Ref. 26. Although the details can be 
improved, the crux of the argument is that the planar density 
of states enters Re- the quasiparticle scatters many times in 
the plane while performing cyclotron motion before hopping 
from plane to plane (measured afe-plane resistivity is of the or- 
der lO/iil-cm as compared il-cm for the c-axis resistivity even 
at optimum doping It is worth noting that oscillations of 
i?c also precisely follows the oscillations of the magnetization 
in overdoped Tl2Ba2Cu06+5.' 
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where (/> = 2'KBa^e/ his a constant. For the matrix B„: 



IV. RESULTS 



(^n)m,m 
(Bn) 

(y-^n^m.m— 1 



-t 



(-1) 



_^/gi(n-i),/, 



(7) 



For the matrix r„, we note that T„ = Bj-^_^j^. 

The 2M Lyapunov exponents, 7^, of 

limAr^oo(7ArT^)^^^", where Tn = 11^=^ Tj' are de- 
fined by the corresponding eigenvalues Ai = e^'. All 
Lyapunov exponents 71 > 72 > . . . > 72M, are computed by 
a procedure given in Ref. 23. The modification here is that 
this matrix is not symplectic. Therefore all 2M eigenvalues 
have to be computed. The remarkable fact, however, is that 
except for a small fraction, consisting of larger eigenvalues, 
the rest do come in pairs (A, 1/A), as for the symplectic 
case, within numerical accuracy. We have no analytical 



There are clues in the experimentsii that disorder is very 
important. For 15 and 16% doping the slow oscillations in 
experiments, of frequency 290 — 280T, are not observed un- 
til the field reaches above SOT, which is much greater than 
Hc2 < IQT. For 17% doping the onset of fast oscillations 
at a frequency of 10, 700T are strikingly not observable until 
the field reaches 60T. The estimated scattering time from the 
Dingle factor at even optimal doping and at AK is quite short. 

For 17% doping corresponding to fj, — — 0.322t and the 
band structure given in Table U a slight change in disorder 
from Vq — 0.7t to Vq = 0.8t makes the difference between 
a clear observation of a peak to simply noise within the field 
sweep between 60 — 62T, as shown in Fig.|2]and Fig.|3] Since 
in this case Wq — Vs = 0, there is little else to blame for 
the disappearance of the oscillations for fields roughly below 
60T. The results are essentially identical for small values of 
Wo, such as Ombt. 
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FIG. 2: (Color online)The main plot shows the Fourier transform of 
the field sweep shown in the inset. The peak is at 10, 695T. The 
inset is a smooth background subtracted Shubnikov-de Haas oscil- 
lations, as calculated from the Landauer formula for 17% doping as 
a function of 1/B. The disorder parameter is Vb — 0.7t.The band 
structure parameters are given in Table |l] 
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FIG. 4: (Color online)The same plot as in Fig. (2] except for 15% 
doping and DDW order. The parameters are given in Table U 
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FIG. 5: (Color online)The same plot as in Fig. [21 except for 16% 
doping and DDW order. The parameters are given in Table U 



FIG. 3: (Color online)The same parameters as in Fig. [2] but Vo — 
Q.dtt. The background subtracted conductance is simply noise to an 
excellent approximation. 



For 15% and 16% dopings we chose Vb to simulate the 
fact that oscillations seem to disappear below SOT. The field 
sweep was between 30 — 60T. The results for DDW order 
are shown in Fig. |4]and Fig. |5] The most remarkable feature 
of these figures is that disorder has completely wiped out the 
large electron pocket leaving the small hole pocket visible. To 
emphasize this point we also plot the results for 15% doping 
but with much smaller disorder Vb = Q.2t; see Fig.|6] Now we 
can see the fragmented remnants of the electron pocket. With 
further lowering of disorder, the full electron pocket becomes 
visible. It is clear that disorder has a significantly stronger ef- 
fect on the electron pockets than on the hole pockets. This, 
as we noted earlier, is largely due to higher density of states 
around the antinodal points, which significantly accentuates 
the effect of disorder'^ We have done parallel calculations 
with SDW order as well. The results are essentially identi- 
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FIG. 6: (Color online) The same plot as in Fig |4l except that 
Vb = 0.2/: instead of 0.8/:. There is now a fragmented electron pocket 
centered around 2100T and the main peak is at ISST. The rest of 
the parameters are given in Table [H 
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FIG. 7: (Color online) The same plot as in Fig. |4]for 15% doping 
but using SDW order. The main peak is at 195T. The rest of the 
parameters are given in Table |l] 
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FIG. 8: (Color online) The same plot as in Fig. [7] except for 16% 
doping and using SDW order. The main peak is at 173T. The rest of 
the parameters are given in Table |l] 



cal. They are shown again for 15 and 16% doping in Fig. [7] 
and Fig. [8] We have kept all parameters fixed, while adjust- 
ing the the SDW gap to achieve as best an approximation to 
experiments as possible. 

It is important to summarize our results in the context of 
experimental observations. First, we were able to show that 
the electron pocket frequencies are strikingly absent because 
of disorder and the slow frequencies corresponding to the hole 
pocket for 15% and 16% doping damp out below about 30 T, 
even though Hc2 is less than 10 T. Similarly, that the high fre- 
quency oscillations at 17% doping do not arise until about 60 
T has a natural explanation in terms of disorder, although in 
this case some magnetic breakdown effect, which was not ex- 
plored, can be expected. This requires both further experimen- 
tal and theoretical investigations. The calculated frequency of 
the high frequency oscillations, 10, 695 T is remarkably close 
to experimental value of 10, 700 ± 400 T. As to the magni- 
tude of the slow oscillations, the calculated values are given 
in Table U which are reasonable in both magnitude and trend 
when compared to experiments. The small discrepancies in 



the magnitude of F are due to our neglect of t" in the band 
structure. This can be, and was, checked by checking the pure 
case, that is, without disorder. 



V. CONCLUSIONS 

In the absence of disorder or thermal broadening, the oscil- 
lation waveforms are never sinusoidal in two dimensions and 
contain many Fourier harmonics. At zero temperature moder- 
ate disorder converts the oscillations to sinusoidal waveform 
with rapidly decreasing amplitudes of the harmonics. Further 
increase of disorder ultimately destroys the amplitudes alto- 
gether. Many experiments exhibit roughly sinusoidal wave- 
form at even ultra low temperatures, implying that disorder 
is important. The remarkably small electronic dispersion in 
the direction perpendicular to the CuO-planes cannot alone 
account for the waveform. 

For NCCO it is no longer a mystery as to why the fre- 
quency corresponding to the larger electron pocket is not ob- 
served. As we have shown, disorder is the culprit. Neither is 
the comparison with ARPES controversial,'"* as in the case of 
YBCO, since there is good evidence of Fermi surface cross- 
ing in the direction (tt, 0) — t- (tt, tt), which is a signature of 
the electron pocket. The crossing along (tt, tt) — > (0, 0) can 
be easily construed as an evidence of a small hole pocket for 
which half of it is made invisible both from the coherence 
factors and disorder effects.''' For electron doped materials, 
such as NCCO and PCCO, it is knowr>ii that the Hall co- 
efficient changes sign around 17% doping and therefore the 
picture of reconnection of the Fermi pockets is entirely plau- 
sible, with some likely magnetic breakdown effects. The real 
question is what is the evidence of SDW or DDW in the rele- 
vant doping range between 15% and 17%. From neutron mea- 
surements we know that there is no long range SDW order 
for doping above 13.4%?^ We cannot rule out field induced 
SDW at about SOT. For DDW, there are no corresponding 
neutron measurements to observe its existence. Given that 
DDW is considerably more hidden '^•^'* from common exper- 
iments, it is more challenging to establish it directly. NMR 
experiments in high fields for suitable nuclei can shed light 
on this question. The unavoidable logical conclusion from the 
QO measurements is that a density wave that breaks transla- 
tional symmetry must be present. We suggest that motivated 
future experiments will be necessary to reach a definitive con- 
clusion. Finally, at the level of mean field theory we have 
been unable to decide between SDW and singlet DDW. At the 
moment the best recourse is to experimentally look for spin 
zeros in the amplitude of quantum oscillations in a tilted mag- 
netic field. A theoretical discussion of this phenomenon that 
can potentially shed light between a triplet order parameter 
(SDW) and a singlet order parameter, the singlet DDW dis- 
cussed here, was provided recently.^^ So far experiments are 
in conflict with each other in YBCO: one group suggests a 
triplet order parameter^" -" and the other a singlet order pa- 
rameter.-^^ 

It is unquestionable that the QO experiments are likely to 
change the widespread views in the field of high tempera- 
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ture superconductivity. Although the measurements in YBCO 
are not fully explained, the measurements in NCCO appear to 
have a clear and simple explanation, as shown here. However, 
given the similarity of the phenomenon in both hole and elec- 
tron doped cuprates, it is Ukely that the quantum oscillations 
have the same origin. 
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